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We establish a fundamental connection between quantum criticality of a many-body sys¬ 
tem, such as Bose-Einstein condensates, and its capacity of information-storage and pro¬ 
cessing. For deriving the effective theory of modes in the vicinity of the quantum critical 
point we develop a new method by mapping a Bose-Einstein condensate of IV-particles onto 
a sigma model with a continuous global (pseudo) symmetry that mixes bosons of different 
momenta. The Bogolyubov modes of the condensate are mapped onto the Goldstone modes 
of the sigma model, which become gapless at the critical point. These gapless Goldstone 
modes are the quantum carriers of information and entropy. Analyzing their effective theory, 
we observe the information-processing properties strikingly similar to the ones predicted by 
the black hole portrait. The energy cost per qubit of information-storage vanishes in the 
large-./V limit and the total information-storage capacity increases with N either exponen¬ 
tially or as a power law. The longevity of information-storage also increases with N, whereas 
the scrambling time in the over-critical regime is controlled by the Lyapunov exponent and 
scales logarithmically with N. This connection reveals that the origin of black hole in¬ 
formation storage lies in the quantum criticality of the graviton Bose-gas, and that much 
simpler systems that can be manufactured in table-top experiments can exhibit very similar 
information-processing dynamics. 
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I. INTRODUCTION 

Perhaps the most mysterious thing about black holes is their capacity to store and process 
information. Already examining the behavior of black hole entropy [I] in the classical limit, one 
arrives at the following puzzle |2j: according to the Bekenstein formula, the black hole entropy 
scales as S ~ where R is the gravitational radius and Lp is the Planck length, which in terms 
of Newton’s coupling Gjy and the Planck constant h is defined as L 2 P = hGpr- This expression 
diverges in the classical limit, if R and Gn are kept fixed. This is because Lp —> 0 in this limit. 
Thus, a quantum physicist, to whom the classical physics is defined as a limit K 0, concludes that 
the entropy (and thus the information-storage capacity) of a classical black hole of a finite mass is 
infinite. On the other hand, a classical observer would seem to have a very different perspective. 

A classical physicist knows that information can only be encoded in the black hole in features 
that are possible to re-arrange. But it is well known that classical black holes are essentially 
featureless. They are uniquely characterized by their mass, charge and the angular momentum. 
Consequently, the classical observer would conclude that black holes should carry zero entropy. 
Thus, it seems that the two observers must disagree on what happens in the h = 0 limit. How can 
infinity and zero be reconciled? 

The above puzzle, along many others, can only be resolved within a microscopic theory, which 
explicitly identifies the quantum degrees of freedom that act as carriers of the black hole informa¬ 
tion. We shall focus on a microscopic theory of this type offered in [3j, in which the black hole 
is described as a composite multi-particle state of soft gravitons at a quantum critical point [2]. 
This system can approximately be described as a Bose-Einstein condensate (or a coherent state) 
of gravitons of characteristic wavelength ~ R and occupation number N R 2 /L 2 P . 

According to this picture, the microscopic carriers of black hole entropy and information are 
collective excitations of the Bogolyubov-type [2]. The remarkable property of these modes is that 
around the critical point they become nearly gapless and decouple from the rest of the modes in 
the N —> 0 limit. Hence, if a system has N gapless Bogolyubov modes, there are expected to exist 
~ e N nearly-degenerate quantum states that can be labeled by the occupation numbers of these 
modes. These states are crowded within a ~ 1 /N energy gap. Correspondingly, the entropy should 
scale as N. 

Such a behavior uncovers a very peculiar nature of black hole information storage and in par¬ 
ticular addresses the above puzzle. In the classical limit N —> oo and entropy is indeed infinite, 
but the energy gap closes to zero. Despite the fact that the number of quantum states located 
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within this gap diverges, the required time for resolving their differences also becomes infinite. This 
time-delay reconciles the views of the two observers. The classical black hole appears as featureless 
for any finite interval of the observation time. 

In the present paper we shall focus our attention on improving the understanding of the relation 
between quantum criticality and physics of information-processing. Since, according to the black 
hole iV-portrait the underlying fundamental ingredient of black hole information-processing is a 
self-sustained quantum criticality of the cold graviton gas, the natural question to ask is: How 
important is the precise nature of the constituent bosons? Of course, for recovering the standard 
gravitational properties of the system, such as the emergence of the Einsteinian metric in a mean 
field description, the spin-2 nature of the constituents is crucial. However, one may suspect that the 
information-processing features that originate from the fact of criticality, may not be so sensitive to 
the precise nature of the constituents. If this is true, then a wide class of critical systems must share 
some of the black hole information-processing abilities. This would be a very important conclusion 
because of the following reasons. First, it would enable us to study black hole properties on much 
simpler systems, both theoretically and experimentally. Secondly, it would open up a possibility 
of implementing the black hole information-processing methods in table-top experiments. 

Indeed, it was noticed already in [2] that the simplest prototype models of attractive Bose- 
Einstein condensates exhibit some striking similarities with the black hole quantum portrait. For 
example, the cost of the information-storage per qubit is by a factor of 1/N cheaper relative to 
the energy-cost exhibited by usual (non-critical) quantum systems. In the latter systems, for the 
same amount of information storage, one typically pays the energy price of the inverse size of the 
system, ~ Also, the degeneracy of states within the 1/N energy gap increases with N. 

Moreover, it was shown that a sharp increase of one-particle entanglement takes place near the 
critical point [5j and the scrambling of information becomes maximally-efficient [6j. The time- 
dependent evolution of the critical condensate uncovers a scaling solution in which the condensate 
is stuck at the critical point throughout the collapse [7|. This is the behavior that one would 
expect if the microscopic foundation of Hawking radiation were through the collapse and quantum 
depletion of the condensate. Q 

In the present paper we would like to take one more step in understanding the universal role of 
quantum criticality for information-processing. To this end, we invent a new method for deriving 

1 Let us remark here that there is interesting evidence from the S'-matrix that the above multi-particle picture 
represents the correct microscopic description of black hole physics. Indeed, the 2 —> N scattering S’-matrix 
element of two trans-Planckian gravitons into N soft ones reproduces an entropy suppression factor e~ N when 
the number N is given by the quantum critical value [3]. Similar conclusions were reached in [5] in a different 
approach. For other aspects and related work, see m 
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the effective low-energy theory of modes near the critical point. As we shall see, the resulting 
effective theory indeed describes nearly-gapless modes with 1/iV-suppressed couplings. This form 
immediately reveals the key features of information storage, such as the low energy cost and the 
long timescale for retrieval as well as the high speed of information scrambling in the over-critical 
unstable phase. 

In our approach we formulate an effective description of quantum criticality in the isomorphic 
language of a pseudo-Goldstone phenomenon in a U[2k + l)-sigma-model. The global U(2k + 1)- 
symmetry of the sigma model corresponds to a continuous unitary transformation among the 
creation and annihilation operators of the bosons with different momenta 0, ±1, ±2,... ± \k\. This 
identification allows us to map the Hamiltonian of the Bose-gas onto a potential of the JJ{2k + In¬ 
signia model. In this way, the problem of diagonalization of the Hamiltonian becomes isomorphic 
to finding the symmetry breaking patterns in the sigma model. The latter is achieved by the 
straightforward minimization of the potential. The resulting pattern determines the spectrum 
of pseudo-Goldstone excitations. The quantum critical points of the Bose-gas correspond to the 
transitions between the different symmetry breaking patterns in the sigma model, where some 
Goldstone modes become gapless. We show that this transition takes place when the particle 
number N and the strength of the attractive interaction a exactly balance each other, aN = 1. 

This approach allows to derive an effective Hamiltonian of the gapless excitations in form of the 
potential energy part of the Goldstone chiral Lagrangian. This form reveals that near the critical 
point the Goldstones represent gapless excitations with 1/iV-suppressed self-interactions, 

.2 


H eff - 


n 


gold 


N 


+ ..., 


( 1 ) 


where n go id is the Goldstone number operator. 

Analyzing the resulting effective theory, we observe a striking connection with anticipated fea¬ 
tures of black hole information processing. Namely, the energy cost per information qubit goes 
to zero as a power of 1/N and correspondingly near the critical point the time for information 
retrieval increases with N. 

The number of states within the collapsed energy gap increases either as a power law or ex¬ 
ponentially with N depending on the nature of the attractive coupling. For a non-derivative self 
interaction the increase follows a power law, whereas for a gravity-type derivatively-coupled case 
the degeneracy is increasing exponentially with N. 

On the other hand, the instability of the system in the over-critical regime implies a positive 
Lyapunov exponent which, as shown in |6j, is the engine for scrambling of information. The log- 
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arithmic scaling of the scrambling time for black holes was originally conjectured in [12] . The 
observation that critical Bose-Einstein systems are fast scramblers provides the microscopic under¬ 
standing of this phenomenon in terms of a Lyapunov exponent and the high density of states near 
the quantum critical point [ 6 |. 

In summary, the effective theory of information at quantum critical point developed in the 
present paper confirms, in accordance with mum, that critical multi-boson systems are indeed 
efficient storers and processors of information. 

Finally, in identifying the information carriers for black hole entropy, we distinguish two dis¬ 
tinct sources of quantum degeneracy. Since degenerate states are labeled by occupation numbers 
of gapless Goldstone modes, we can create degenerate states either by increasing the occupation 
numbers of a small subset of Goldstone modes, or by diversifying the population among many 
different Goldstone modes, while keeping individual occupation numbers small. The crucial differ¬ 
ence among the two types of states is that the information stored in un-diversified large occupation 
numbers can in principle be resolved in interference experiments even in the classical limit. Irre¬ 
spectively whether contributions from such states must be counted into black hole entropy, their 
existence gives a very interesting possibility of detecting the quantum sub-structure in macroscopic 
experiments. 


II. UNDERSTANDING THE CRITICAL POINT IN TERMS OF A 
PSEUDO-GOLDSTONE PHENOMENON 

We start with the system that was used in [2] as the simplest prototype model for exploring 
the connection between black hole information processing and quantum criticality. It is a gas of 
cold bosons with a delta-function-type attractive interaction. Such systems are well known and 
have in particular been studied in the context of cold atoms [10] . We consider a Hamiltonian of 
the following form, 

n = [ d d x^+^^ 

J 2 m 

where t/j = -Lg '* 1 y = R d j s the d-dimensional volume and k is the d-dimensional 

wave-number vector, at, a^ are creation and annihilation operators of bosons of momentum k. 
They satisfy the usual commutation relation [a^, at,] = . g is the coupling. Rescaling the 

Hamiltonian, we can write % = 2R? m ^• Introducing a notation a = (yji) ^ 7 ^ and taking d = 1 


^ — gh / d d x p + . 


( 2 ) 





we arrive at the following Hamiltonian, 
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H = Y k2 4 °* - \ Y 4 , 42 ^ 3^4 • ( 3 ) 

k ki+k2—ks—k4=0 

We restrict ourselves to three levels, k = 0, ±1. Legitimacy of this approximation shall be justified 
later. 

Let us define a triplet operator a* = (a_i,ao,ai). Correspondingly, we define the number 
operators n* = aj a t . We treat the operator a* as a triplet under a global symmetry group 517(3), 

ai a' = 17/, (4) 

where 17/ is a unitary transformation matrix which keeps invariant the total number operator, 

n = Y n i- ( 5 ) 

i 

Thus, in any state with non-zero particle occupation number, 


<n) = N, 


( 6 ) 


the SU (3)-symmetry is spontaneously broken and the order parameter of breaking is N. This fact 
allows us to map the dynamics of the quantum phase transition in a Bose-Einstein condensate 
described by the Hamiltonian ([3]) onto the pseudo-Goldstone phenomenon in an 517(3) sigma 
model. The Hamiltonian in terms of the triplet ai can be written as 

l 

ni - - (a[al_ia 0 a 0 + a' 0 a' 0 aia-i) + t 2, "1 + H SU(3) . (7) 

*=-i,i 4 »=—l 

where, Hgu ( 3 ) is the 517(3)-invariant part of the Hamiltonian and has the following form, 

a o ot 

H S u( 3) = - ~n 2 + - n + n(n - N). ( 8 ) 


Here, /x is the Lagrange multiplier that fixes the total particle number. Let us first investigate 
the ground-state of the 517(3)-symmetric Hamiltonian. Minimization with respect to rii gives the 
following equations, 


dH S u( 3 ) 
drii 


a 


— an + — + 11 = 0 
4 


(9) 


and 


dH S u(z) 

dfi 


= n — N = 0 . 


( 10 ) 
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The vacuum is achieved at n = N and fi = a(N — 1/4). The SU (3)-symmetry is spontaneously 
broken down to SU{2) and there is a doublet of Nambu-Goldstone boson^j For example, choosing 
the vacuum at n\ = n_i = 0, 713 = N and ignoring 0(l)-corrections to the zero mode occupation, 
the Goldstone doublet is ( 01 , a_i). The existence of the Goldstone boson is the manifestation of the 
indifference of the system with respect to redistributing occupation numbers among the different 
levels. 

Now, the addition of the first three terms in ([7]) breaks the SU (3) symmetry explicitly and lifts 
the vacuum degeneracy. Let us investigate the effect of these terms on the state n* = (0,1V, 0). 
Notice, this VEV continues to be the extremum of the Hamiltonian even after the addition of the 
explicit-breaking terms, but the value of /j is now shifted to /r = §(-/V — 1/2). The would-be 
Goldstone modes now acquire a non-trivial mass-matrix, 


1 - - N 

-L 2 ’ 

- -N 
2 ’ 



( 11 ) 


This matrix is diagonalized by states au g ht = ^(°i + a -i) and a heavy = ^j( a i — a -i)) with 
eigenvalues 1 and 1 — aN respectively. In this basis the Hamiltonian becomes, 


1 , 0 \ 

( 12 ) 

0,1 — aN .J 

Thus, a would be Goldstone doublet is now split into a light mode and its heavy mode partner. 
Since the explicit breaking terms leave unbroken the U( 1) subgroup of 517(3), corresponding to the 
generator Q = diag( 1,0, —1), the mass eigenstates a ugh t, and aheavy are also U (l)-charge eigenstates. 
They carry charges equal to +1 and —1 respectively. In the usual language, the charge Q is a 
momentum operator. 

Let us at this point highlight the relation between the Goldstone modes au g ht and aheavy and the 
usual Bogolyubov eigenmodes of the mass matrix ©• The latter are obtained by diagonalizing 
the mass matrix through a canonical transformation; in our setup, they are given by 


6-hi = itia±i — uia- Fl 


(13) 


with 


1 + y/1 — aN 1 — y/1 — aN 

2(1 - aN) 1 / 4 ’ Vl ~ 2(1 - aN") 1 / 4 ' 


(14) 


2 When counting the number of physical Goldstone modes, we have to take into account that only different rays in 
the Hilbert space are independent. A breaking SU(N) —> SU(N — 1) can thus be understood as the breaking of 
SO(N ) —» SO(N — 1) on the Hilbert space and gives rise to N — 1 Goldstone modes. 
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In this description, the Hamiltonian has two degenerate eigenvalues 

ei_i = Vl- aN . (15) 

The modes ai; g ht and aheavy relate to the Bogolyubov modes through 

“■«“ = ^ u ' +n) ( b ' + i, h) = ^ (1 ' aN)1/i (‘i+ 6 b) . (“) 

«hei»T = ^(“1 - Hi) (&1 - ill) = —-- (l>l - til) - (I?) 

The Goldstone mode thus encodes a light direction in configuration space, reachable through 
occupying both Bogolyubov modes simultaneously. Note that a state with Bogolyubov occupation 
number of 0(1) corresponds to a Goldstone configuration with «-j ig } lt /hight ~ 1 /\/1 — aN, which is 
responsible for the y/1 — alV-difference in the eigenvalues. 

The heavy partner, on the other hand, corresponds to a heavy multi-particle direction in 
configuration space, reachable through populating either one of the Bogolyubov modes with 
0 { 1 /vT — adV) particles. In turn, this explains the alV-independence of the second eigenvalue of 

& 

For 1 — aN < 0 one of the pseudo-Goldstone bosons becomes tachyonic and the vacuum is 
destabilized. The physics of this instability is that for aN > 1 the SU( 2 ) preserving vacuum is 
no longer energetically favorable and the system flows towards the new ground-state. This new 
ground-state can be found by minimizing the full Hamiltonian, 0. In order to minimize this 
Hamiltonian, let us set a* = s/nte™*. The only phase-dependent term in the Hamiltonian is the 
second term in 0, which takes the form 

- ^ (ala^aoao + ajajaia_i) = -a ( A /nin_in 0 )cos( 6 'i + 6 Li - 2 8 0 ). (18) 

Since there is no other conflicting phase-dependent term in the energy, in the minimum we will 
have cos(0i + 0_i — 2#o) = 1- We can thus set without any loss of generality, 

9 \ + 0_i — 28 q = 0 . (19) 


Another simplifying observation is that in the extremum n\ and n_i must be equal. This can be 
seen by extremizing the Hamiltonian with respect to n\ and n_i, 


dH a In -1 a a 

d^ = 1 ~ 2VV n ° + 2 ni ~ an+ 4 + M = ° 


( 20 ) 


dH 

dn-i 


= 1 - 


n\ a a 

-no + 7 T n -i ~ an + — + /x = 0 

n_ i 2 4 


a 

2 


( 21 ) 
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which show that n\ and n_i can be non-zero only together. Moreover they must be equal, since 


they satisfy the same equation, with only one positive root. This can be seen by multiplying (20) 


and (21) by n\ and n_l respectively. The resulting quadratic equations are identical, 


dH a 9 a a _ 

ni -— = -n\ + ni(l - aN + -+//)- -y'n.inino = 0 

OTi\ Z 4 z 


( 22 ) 


dH a 9 a a _ 

n-1- - = — n_, + n_i(l - aN +-+//) - -^/n_inin 0 = 0 

uTi—i Z 4 Z 


(23) 


and have only one positive root. Thus, without any loss of generality we can minimize the Hamil¬ 
tonian for the configuration n* = (x, N — 2x,x), which gives 
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a 


H = -ax 2 + 2x(l - aN) - jN (N - 1) . 


(24) 


Since x is positive definite, for aN < 1, the minimum is achieved at x = 0. For aN > 1, the 
minimum is at 


X = Ya^ aN ~ ^ 


(25) 


Not surprisingly, this corresponds to the first Fourier modes of the bright soliton solution to the 
Gross-Pitaevskii equation in the overcritical regime jTOlj . The energy of the ground-state is given 
by 


Hmin = - — {aN - If - —N{N - 1), 
7a 4 


(26) 


where the first term only exists for aN > 1. 

We can now understand the critical phenomenon in terms of this Goldstone-mode. The large 
occupation number serves as an order parameter for a spontaneous breaking of a global SU(3)~ 
symmetry. This symmetry corresponds to redistribution of the particle occupation numbers be¬ 
tween the different momentum states, without changing their total number. This symmetry is 
only approximate, and normally the would-be Goldstone mode has the mass of the order of the 
first momentum level (that is order one in our units). However, the phase transition corresponds 
to the point where the particle number distribution changes. This necessarily implies that the 
corresponding pseudo-Goldstone mode must become massless at this point. 

Notice that since x parameterizes the occupation number of (pseudo-)Goldstones, x = n go id = 


a \oid a 9 oidi the Hamiltonian (24) essentially represents the effective action for this mode, 


a 


HQold, flgold) Qgold, T ^goldJ^gold < ^ -0 • 


(27) 
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FIG. 1. (a) Ground state for aN = 1.5 at various values of N, using the parametrization | GS) = 

J2x=o a x\%, N — 2x, x ). a x are normalized according to J2x=o a x = N, since as N —> oo, J^ 2 a(Ny) 2 dy = 1 
and thus it’s possible to compare them at different N. As N increases, the distribution becomes increasingly 
more peaked at x m . (b) x m for N = 800 particles compared with the analytic prediction x m = 1). 


Thus, the mass of the pseudo-Goldstone is given by m 2 gold = 2(1 — aN) and the self-coupling is 
given by a go id = a\. At the critical point, the Goldstone mass term vanishes and the theory is 
described by a gapless mode with a self-interaction strength given by a go id = 7/2 N, 


17 1 

Hooid = (n g oid) 2 ^2 “ 1 N ■ 


(28) 


Of course, the phenomena uncovered in the Goldstone formulation are in one-to-one correspon¬ 
dence to those that are seen in a mean-field and Bogolyubov description of the model. In fact, 


as hinted before, the minima of the Hamiltonian Eq.(24) correspond to the k = —1,0,1 modes 


of the solutions to the mean-field Gross-Pitaevskii equation: for aN < 1 we found the homoge¬ 
nous condensate, while the solution for aN > 1 corresponds to a localized bright soliton. In this 
context, one should also keep in mind that the solution (x, N — 2x, x) only agrees with the exact 
quantum mechanical ground state in the limit N —> oo. At any finite N, the ground state is a 


smeared distribution centered around (x,N — 2x,x) with x given by (25), as can be seen in Fig. 
([I]). In particular, the ground state at the critical point is characterized by strong entanglement 
and is therefore not well described by a mean field. This, however, will not alter our conclusion on 
state-counting and information processing to be performed in the following sections. 

We reemphasize that in the Goldstone language, the critical point amounts to destabilization of 
the S'?7(2)-invariant vacuum. With the pseudo-Goldstone method we have traded the diagonaliza- 
tion of the Hamiltonian with the minimization procedure. Substituting the small deformations of 
the order parameter by the Goldstone mode allowed us to derive the effective action for the latter. 
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This is similar to deriving an effective Hamiltonian of a phonon field in a background external 
magnetic field, which breaks the rotational symmetry of a ferromagnet. The analogous role in our 
treatment is assumed by the terms in the Hamiltonian that explicitly break SU (3)-symmetry. [^] 
We conclude this section by pointing out that one can continue an analytic treatment of the 
system into the solitonic regime by performing an x-dependent canonical rotation on the creation 
and annihilation operators a,. We define the condensate mode 


t-cond \l AT 1 T \l 1 n + \j AT ®1 j 


and the corresponding orthogonal modes 


c±i = 


±1 + 1 1 - 


2x 


n 


a_ i — 2 


x 1 

JV° 0+ 2 


+ 1 + 1/1- 


2x 


n 


a~ i 


(29) 


(30) 


The resultant Hamiltonian is now always minimized by the configuration (c_i, c con d, ci) = (0, N, 0). 
Expanding around this vacuum leads to a quadratic Hamiltonian which can be written in matrix 
form as 


H‘2 = (c^cl^CljC-l 


( m\ 772-2 7723 7724^ 


7712 7771 7/74 7773 
7773 7774 777i m 2 
\m^ 7773 7772 m \/ 


( Cl \ 


C-l 

c\ 

V c -1 ) 


(31) 


with 


/ 772i\ 

772 2 

772 3 

\t72 4 7 


(l-^ + x(8a-?±]^)\ 
(—1 + 3aN — 7ax) 

ax (o 7 x \ 

2 V' 3 

“ f 4x - A' - 




2V 


(32) 


Canonical diagonalization again leads to two distinct eigenmodes; for aN < 1, their energy is 
e = yd — oiN. For aN > 1, one of them remains exactly gapless, while the other has a frequency 
e = \/l — aN y/|(4 + 3aN). The former is the mode that corresponds to translations of the soliton, 
whereas the other is the Bogolyubov mode that is only light in the vicinity of the critical point. 
In the Goldstone language, the Goldstone mode aught splits into two parts, essentially given by its 
real and imaginary parts. One remains gapless and generates translations of the localized ground 


Notice also a curious analogy with the sigma-model of large-IV-color QCD. In this case our N would be mapped on 
the number of colors, whereas the levels k = ±, 0 to quark flavors. The pseudo-Goldstone boson then is analogous 
to pion, which also has 1/N suppressed self-coupling. 
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state; the other is gapless only at the critical point and obtains a frequency e = (1 — aN) 
for aN > 1. 

The appearance of the light modes can of course again be related to the breaking of SU( 3) 
generators. Keeping in mind that in the rotated description the symmetry transformations are given 
by d = SUS~ 1 c, where S is the symplectic matrix that generates the canonical transformation, 
we immediately see that the translation generator Q = diag( 1,0,—1), which commutes exactly 
with the Hamiltonian for all x , does no longer annihilate the groundstate. The generator that 
redistributes particle numbers, on the other hand, only commutes with H at the critical point. 


III. INFORMATION PROCESSING 


The above form (27) of the effective Hamiltonian displays the role of quantum criticality for 


information storage and processing. The quantum information in the above system is encoded 
in the state of the Goldstone mode. The remarkable thing about it is the low energy cost of 
information-qubit-storage, which is suppressed by powers of 1/N relative to the inverse size of the 
system. This phenomenon is a manifestation of quantum criticality. 

In the time evolution of the Goldstone state we can distinguish two sources. One is the interac¬ 


tion governed by a quartic self-interaction Hamiltonian (28). The rate of this process is suppressed 


by powers of 1/N, and correspondingly the time-scale of evolution is very long. 

The second source of time evolution is the Goldstone mass term that parametrizes the depar¬ 
ture from quantum criticality. In the overcritical regime, this mass is imaginary and results in 
an exponential growth of the Goldstone occupation number. This instability is described by a 
Lyapunov exponent that, as shown in [6], leads to the generation of one-particle entanglement and 
potentially to the scrambling of information. 


IV. NUMERICAL RESULTS OF STATE EVOLUTION 


An exact numerical diagonalization of the Hamiltonian ([3]) provides a complementary analysis 
which is valid also at the critical point. Using the same technique as in |i5], we can verify the above 
results by comparing the the first Bogolyubov state as well as a lowly occupied Goldstone state to 
the exact eigenstates of the Hamiltonian ©• Several quantities are instructive. 

In Fig{2] we plot the expectation value of the energy in a state with Goldstone occupation 


“light“light) = 1 as a function of aN. The solid line corresponds to the analytic result (27). The 
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FIG. 2. Expectation value of energy in state with Goldstone occupation \ <4 ight aught y = 1 

as a function of aN for N = 500. 


D„(t) 


v 




FIG. 3. (a) Normalized exact time evolution of the time dependent part of the first Bogolyubov state for 
N = 500 and aN = 0.004 (black, solid), 0.596 (red, dashed) and 0.976 (blue, dotted), (b) Frequency of 
oscillation as a function of aN at various values of N. The solid line corresponds to the analytic result 
2 a/1 - aN 


increase of energy around the phase transition is suppressed by 1/N and thus not visible. 


In Figj^i we plot the exact time evolution of the Bogolyubov state |1#) = F^I-JOb), where |0 b) 
is the Bogolyubov ground state. The results are obtained for fixed particle number N = 500, while 
the effective coupling aN is varied. We observe an decrease of the frequency as aN approaches 
the critical point. To illustrate this point better, we plot the frequency of oscillations versus aN 
for N = 100, 300 and 500 in Fig [3 ]d. We observe the square root behavior (15), with the frequency 
scaling with aN as /I — aN + 0(1 /IV). 
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V. OCCUPATION AND STABILITY OF HIGHER MODES 


We can now verify that our starting assumption about the stability and a negligible occupation 
number of higher momentum modes is valid as long as < 1. As we shall see, the occupation 
numbers of modes with momentum number k are suppressed by a factor of (x/A r )l fc l“ 1 . This is a 
manifestation of the fact that the global SU(2\k\ + l)-symmetry, corresponding to redistribution 
of occupation numbers among the fist k levels, is explicitly broken by the momentum term in the 
Hamiltonian <©■ We can think of this breaking as the decoupling limit of a spontaneous breaking 
by a spurion order parameter, £, in the adjoint representation of the SU(2\k\ + 1) group, with 
the expectation value £ = diag(—|fc|,..., —1,0,1,..., |/c|)j^] Obviously, the kinetic term in the 
Hamiltonian ([ 3 ]), can then simply be written as an SU(2\k\ + l)-invariant product, a 1 

The spurion £^£ can be thought of as an expectation value of the number operator of a ’’hidden” 
Bose-gas of some (jj-particles, interacting with a* by an infinitely weak coupling n, s.t. £ + £ = 
k'Y2((t In that case, the momentum operator is replaced by a peculiar SU(2\k\ + 1) invariant 
interaction term and both occupation numbers spontaneously break one and the same SU(2\k\ + 1)- 
symmetry. The Goldstones bosons are superpositions of a-s and a-s. The orthogonal combinations 
are pseudo-Goldstones. 

Now, let us take the limit k —> 0, while simultaneously taking the occupation numbers of 
(7-particles to infinity in such a way that the products n(a\a,) = i 2 are kept fixed for all i = 
— 1,1,.... We arrive at the situation in which in the a-particle sector, the entire information about 
the breaking is summed up in the expectation value of the spurion £. Correspondingly, all the 
Goldstones in the a-sector become pseudo-Goldstones. With increasing momentum level \k\, the 
masses of pseudo-Goldstones increase whereas their occupation numbers rapidly diminish. 

In order to see this, let us compute the occupation numbers of the k = ±2 modes. The 
corresponding annihilation operators for k = +2 and k = —2 are 02 and a_ 2 , respectively. The 
bilinear mass matrix of these modes expanded about the state n* = (0,x,N — 2x,x,0), has the 
form 


4 — aN , -fJV 


-flV, 


(33) 

4 - aN. 

This has two eigenvalues, 4 — ^aN and 4 — | aN. Both eigenstates are stable as long as aN < 8/3. 

However, the modes ?r_ 2,2 nevertheless can get populated due to the mixing with the lower 
levels. The occupation number of the lower modes acts as a source for the higher ones. The 

4 This expectation value leaves invariant the [/(l)-subgroup with the generator Q which is proportional to S. As in 
the case of SU(3), this charge represents a momentum operator. 
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FIG. 4. Contribution of the |fc| = 2 modes to the lowest lying eigenstates for IV = 10 (black) and N = 20 
(red). The dots correspond to the exact numerical results, while the solid line depicts the analytical behavior 

(1351- 


occupation numbers of the higher modes can be easily obtained by minimizing the bilinear Harnl- 
tonian including the effective source terms for a_2,2- This source is generated from the interaction 
term after plugging the expectation values of 71 - 1 , 0 , 1 . The phases of the expectation values can 
again safely be set to zero, because all the source terms are negative. The corresponding effective 
Hamiltonian for <22 and a _2 has the form, 

H 2 -2 = ^ ((8 — 3 aN){a\ + a_ 2 )(a 2 + aL 2 ) + (8 —alV)(a^ - a_ 2 )(a 2 - at 2 )) - 

- -axy 7 (N - 2x)(a 2 + at 2 ) + h.c. (34) 

Minimizing this with respect to a_2,2 we see that 

n_2,2 ~x 2 /iV. ( 35 ) 


Similarly, we can derive the occupation numbers of higher momentum modes. In general, modes 
with momentum \k\ have occupation numbers ~ x(x/N)^^ 1 , which rapidly approaches zero for 
x/N <C 1. Thus in this regime, our assumption is well justified. In particular, this regime covers 
the neighborhood of the phase transition, around which x is small. 

This behavior can also be confirmed by a numerical study of the system. To this end, we 
consider the Hamiltonian ([ 3 ]) with a momentum cut-off at \k\ = 2. An instructive quantity is the 
occupation number n_ 2,2 of the \k\ = 2 modes for low lying states around the phase transition. As 
can be seen in Fig. [4j there is no appreciable contribution of the \k\ = 2 modes to the low lying 


states at least until aN ~ 2. We can moreover confirm the quadratic behavior (35). 
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VI. DERIVATIVELY-COUPLED CASE 


The model considered in the previous sections admits two solutions that become degenerate for 
IV —y oo at the critical point. In d dimensions, this number can be increased by a factor of d; this 
corresponds to the number of gapless Bogolyubov modes. For finite A, we may define an entropy 
for a given A-particle system at the critical point by counting all states with < 1/A. Given 

the gap of the light modes, A E ~ 1/A, we obtain J2k ( d+ fc _1 ) ~ A d states within the accountable 
range of energies. This implies an entropy S ~ dlog A. 

We may attempt to reproduce an entropy that scales with A like the Bekenstein-Hawking 
entropy of black holes, S ~ A, by increasing the number of pseudo-Goldstone modes at the critical 
point. This can be achieved, for example, by taking the coupling to be momentum-dependent. In 
order to see this, consider a Hamiltonian, 

H = j d d x'ip + 2 ~ L d+1 h J d d x (^p + 'Vip + ^j , (36) 

where i/j = a*,, V = R d is the d-dimensional volume and k is the d-dimensional wave- 

number vector. L* is a fundamental length and sets the cutoff of the effective theory. Rescaling 
the Hamiltonian, we can write T-L = „ f H , where 

H = Yk 2 a\a k ~ a o y (k 2 k^j ■ (37) 

fc /ci+&2 — &3 — &4=0 

( \ o d 

and ao = ( y r ) Let us again find the effective bilinear Hamiltonian for k p 0 modes, about 

the point no = A and n^o = 0- We obtain 


H = y k 2 (1 - ao A) ata fc 


(38) 




Thus, at ao = 1/A, all modes are critical and there is the same number of massless pseudo- 
Goldstones as the number of momentum modes. 

This peculiar behavior can be equivalently understood in a mean field analysis. Minimizing Eq. 


(36) under the constraint of fixed total particle number yields the Gross-Pitaevskii equation 

L d+1 


A ( 

- ib + 

2 nr 


-'i/i'l’A'0 2 = /a/ j 


(39) 


where /i is again the Lagrange multiplier fixing the total particle number f dVi/j^iJi = A. Using 


the ansatz Un = \ ^e lkx , Eq. (|39|) becomes 


m 


— - - «oA = fi. 


(40) 
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Hence plane waves of arbitrary wavenumber solve the GP equation (39). By inserting ipQ into the 


Hamiltonian (36), we read off the energy of the plane wave solutions 

P 

E k = — (1 - a 0 N) . 

2m 


(41) 


At the critical point, all plane wave solutions become degenerate; the corresponding modes are 
the Bogolyubov modes. 

In order to give a lower bound on the entropy of this system for a given N, we have to keep 
in mind that for the validity of the non-relativistic treatment we should only include modes with 
\k\ < kmax^mR/fo- Their number is N max ~ (k ma x) d - This number sets the number of legitimate 
pseudo-Goldstone modes. However, not all the massless-pseudo-Goldstone modes contribute the 
same weight into the entropy. This is because their self-couplings have different strengths. Hence 
exciting different pseudo-Goldstone species will contribute into the Hamiltonian differently. In 
order to identify which species give the maximal contribution into the entropy of the system, let 
us consider the effective Hamiltonian for the pseudo-Goldstone modes. It has the following form 
(we ignore numerical factors OF order one), 

Hcoid = ^2 \k\ 2 ( (rigoidik)') a 0 + n gM (k ) (1 - a 0 N ) J + cross-couplings, (42) 
0 

where n go id(k) = a go id(k)^a go id(k) is the occupation number of pseudo-Goldstone of a given mo¬ 
mentum number k. At the critical point, the mass terms vanish and the Hamiltonian is given by 
the quartic couplings, 

\ 2 


Hcoid — E \k\ 2 ao (n go id{k2j + cross-couplings. 


(43) 


fc^o 


We can obtain a lower bound for the number of states in the allowed range by looking for the \k\ 
that provides the maximal contribution to the partition sum. We need to take into account two 
competing effects: The larger | k \, the more modes, and thus states, are supplied. On the other hand, 
the energy cost of a given mode grows with |fc|, thereby limiting the maximum allowed occupation. 
Moreover, we pay a penalty for occupying modes more than once as seen from the quadratic 


dependence of (43) on n go id{k). A conservative lower bound is thus obtained by summing over all 
states in a shell around the limiting momentum ki which allows for all Goldstone modes of given 
wavenumber to be occupied exactly once. Lower momenta will yield subdominant contributions due 
to the reduced number of modes and the related penalty for higher occupation; higher momentum 
modes, on the other hand, are too costly to excite. We obtain the contribution, 


SkHcoid 


OL 0 


(44) 
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where we ignored a possible cancellation from the cross-couplings. The factor |/q| rf comes from the 
total number of modes. The maximal contribution to the entropy will come from the species with 
largest \k\, subject to the condition < 1- Taking into the account that «o = 1/-/V, and 

assuming L*/h = m _1 , we obtain 


\ki\ = N d + 2 , n modes = N d + 2 


(45) 


for the limiting momentum and the corresponding number of modes. We get the following contri¬ 
bution to the entropy from the pseudo-Goldstone species, 


IV* 


log y. ( • I ~ n d + 2 


k 


(46) 


Note that we have also evaluated the full partition sum numerically and thereby verified that the 
wavenumbers with |fc| ~ \ki\ provide the dominant contributions. 

If we take into the account the effect of cross-couplings, ignoring possible ’’flat directions” on 
which the different contributions cancel, we get a more conservative lower bound on the entropy. 
The number of cross-couplings for a given momentum level \ k\ scales as the number of corresponding 


pseudo-Goldstones squared, ~ ( kd) 2d . Then, equation (44) changes to 


6 k H Gold ~ \k\ 2d (\k\ 2 a 0 ) 


(47) 


This lowers the allowed number of simultaneously-excitable Goldstones to 


IV* = IV 2d+2 . 


(48) 


Independently, we observe the general message that the derivative self-coupling of bosons leads 
to a dramatic increase of pseudo-Goldstone modes. 


VII. THE CARRIERS OF BLACK HOLE ENTROPY 

We finally wish to address the following question. Which of the large number of degenerate 
states are good candidates for carrying the Bekenstein-type entropy of a black hole? In order to 
answer this question, we shall distinguish two categories of states. 

1) The first category represents states that cannot be resolved semi-classically in any macro¬ 
scopic measurement. That is, the quantum information stored in such states becomes unreadable 
in the infinite -N limit. We shall refer to this category of states as type- A. 
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2) The second category are states that can be resolved in some macroscopic interference exper¬ 
iments. That is, the quantum information stored in such states can be read out even for N = oo. 
We shall correspondingly refer to these states as type-B. 

What distinguishes these two categories of states microscopically? In our picture the states of 
both categories can be labeled by the occupation numbers of some nearly-gapless quantum degrees 
of freedom. As discussed above, these information-carriers can be described as Bogolyubov and/or 
Goldstone degrees of freedom. What distinguishes the two category of states is the scaling behavior 
of information-carrier occupation numbers in large N limit. 

The type-A category of states refers to those in which the relative occupation number of gapless 
degrees of freedom vanishes in the large -N limit. That is, none of these information-carrier degrees 
of freedom are macroscopically occupied. Correspondingly, the type-A states cannot be resolved 
in any macroscopic measurements. 

For type-B states, some of the Goldstone modes can be macroscopically occupied with an 
occupation number that scales as a non-vanishing fraction of N in the N —> oo limit. Such states 
can be resolved in macroscopic interference experiment^} 


This discussion is more relevant for the derivatively-coupled model (36), since it exhibits a large 
diversity of Goldstone modes. However, to keep the discussion as simple as possible, we shall here 
consider only the model ([2]). Therein, let us now exemplify which states can be attributed to the 
type-A and type-B categories. 

As discussed at the end of section [Hj at the critical point, two gapless modes emerge. The 
appearance of these states is a direct consequence of quantum criticality. They come from those 
components of the SU (2)-doublet pseudo-Goldstone that are not affected by the explicit breaking. 
One of them is only an approximate Goldstone, corresponding to one linear superposition of the off- 
diagonal generators from the quotient SU (3) / SU (2) x U( 1). This Goldstone is massless only near 
the critical point. In contrast, the other degree of freedom is an exact Goldstone corresponding to 
the spontaneously broken [/(l)-generator Q = diag( 1, 0, —1) of 577(3). This Goldstone is massless 
everywhere at and beyond the critical point. 

A counting of states in the Goldstone as well as in the Bogolyubov language suggests that the 
number of independent states below the 1/N energy gap near the critical point scales at least as 
N 1 / 4 . However, the total number of almost orthogonal states with expectation value of energy be¬ 
low the same gap becomes larger if we also include coherent states in the counting (see Appendix). 

5 In Appendix [a] we connect to the classical limit by performing the counting of states in a coherent state basis. In 
this basis, one can immediately write down the matrix elements that describe macroscopic interference experiments 
capable of resolving differences between the different type-B states. 
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This is because the coherent states are neither energy nor number eigenstates and involve superpo¬ 
sitions of arbitrary energetic modes with arbitrarily large occupation numbers. Correspondingly, 
they explore a much bigger fraction of the Hilbert space than the states constructed out of a finite 
number of energy or number eigenstates. 

Only the states with vanishing fractions of Goldstone (or Bogolyubov) occupation numbers 
belong to the type-H category. On the other hand, states with macroscopic occupation numbers 
of Goldstones, such as coherent states, belong to the type-B category. The latter states can be 
resolved in macroscopic measurements, even in N = oo limit. This is a manifestation of the fact 
that the coherent states are “classical”. 

Given the fact that interference experiments with black holes have never been performed, both 
of these categories are very interesting, since they both reveal the internal microscopic structure 
of black holes. However, the question is whether both of the types of states should be counted as 
carriers of black hole Bekenstein entropy. 

Should Bekenstein entropy originate exclusively from the quantum states with small-occupation 
numbers of many species of gapless Bogolyubov/Goldstone degrees of freedom, or should the states 
with large occupation numbers of few species also count? 

For now, we shall keep both options open. 


VIII. CONCLUSIONS 

The purpose of this paper was to emphasize the fundamental role of quantum criticality for 
information-processing in many-body systems and to establish this phenomenon as a viable candi¬ 
date for the underlying mechanism for black hole information-processing, along the lines of El Eli. 
We showed that the simplest multi-particle systems with critical behavior capture qualitative fea¬ 
tures of information-processing that is expected for black holes, such as the low energy cost of 
information storage, large degeneracy of states and scrambling of information. 

In order to visualize the nature of the critical phase transition in the language of spontaneous 
symmetry breaking, we developed a new description in which we mapped the quantum phase 
transition in attractive Bose-gas on a Goldstone phenomenon in a sigma model. The two systems 
represent two realizations of one and same unitary symmetry that rotates different momentum 
modes into each other. This symmetry is broken both spontaneously as well as explicitly, but 
the explicit breaking vanishes, up to 1/N effects, at the critical point, resulting onto the gapless 
pseudo-Goldstone modes. 
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This mapping allows us to establish a correspondence between the quantum criticality in Bose- 
Einstein systems and the change of symmetry breaking pattern in a sigma-model. Correspond¬ 
ingly, we map the entropy-carrying (nearly)gapless Bogolyubov modes of the Bose-gas on the 
(nearly)gapless Goldstone modes of the sigma model. This mapping also allows to formulate, up 
to 1/fV-corrections, an effective theory of nearly gapless modes in the regime in which occupation 
numbers are large. 

Our findings, within the validity-domain of the description, confirm the results of previous 
studies mm® and shed light on the criticality phenomenon from a novel angle. Our studies 
indicate that the key information-processing properties of black holes are shared by a larger class 
of the critical systems, including the ones that can be designed in table-top labs. This opens up 
a possibility of “borrowing” black hole information processing abilities for implementing them in 
the laboratory systems. 

Since our findings suggest that critical instability is the key for efficient information processing, 
it would be interesting to find out, by generalizing our approach, whether other unstable systems, 
for example, the ones with parametric resonance instabilities [13] and large classical statistical 
fluctuations, studied in m, also exhibit some analogous properties of information scrambling and 
processing. 
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Appendix A: Counting in Coherent State Language 


In order to highlight the difference between the above two types of states, it is useful to perform 
the state-counting in a coherent-state basis. We shall do this on the simplest example of system 
Q. Generalization to the derivatively-coupled case is straightforward. 

The reason for choosing a coherent-state basis is due to its usefulness for establishing the 
contact with the classical limit. Since the characteristic property of typ e-B states is precisely 
their distinguishability in the classical limit, the coherent state picture is more convenient for their 
description. In particular, in this basis we can immediately write down the matrix elements that 
describe macroscopic interference experiments capable of resolving differences between the different 
type-F> states. 

Let us evaluate the expectation value of the three-level hamiltonian ([3]) over a set of coherent 
states with the fixed total average number N. These states can be labeled by parameter x and the 
three phases 0o,#+i,6Li respectively. We shall denote them by |a;, do, 9 + \, 9-\). Each such state 
represents a tensor product \x, #o, #+i, #-i) = rL-=o ±1 ®\Nk)coh of three coherent states, one per 
each momentum level k. These have the following generic form, 


n k 

N k “ N, 2 e lJlk9k 
I N k ) coh = e~* Y, /z-r K> 
o 


yMk'- 


(Al) 


where k = 0,-1,+1. The mean occupation numbers satisfy Nq = N — 2x and iV_i = N\ = x. 
The integers = 0,1, ...oo label the Fock space occupation numbers for each of the three levels. 
Taking the expectation value of the Hamiltonian and imposing the minimum constraint on the 


phases (19), we arrive (not surprisingly) at the effective Hamiltonian (24). 


Obviously, we should only count the coherent states that form an almost-orthogonal set. The 
scalar product of the two distinct coherent states | tp) = |x, $+i> $-i) and | ip') = \x', 9' 0 ,9' +1 , 


is, 


\{^')\ 2 = e -S fc =o,±i((v^-v^) 2 + 2 \/j^(i-coB(flfc-^))) _ 


(A2) 


Thus, the suppression of the matrix element between the two coherent states comes from the two 
sources: The difference in occupation numbers and difference in phases. 

We shall only be interested in the set of states over which the expectation value of the Hamilto¬ 
nian stays within ~ 1/N gap. For such states the value of x is bounded from above by x max ~ V~N. 
Thus, without taking into the account the variation of phases, we have maximum ~ A 1 / 4 coher¬ 
ent states that are orthogonal due to difference in the occupation number x. In addition, there 
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is a degeneracy due to variation of the phases. For a fixed x, we can perform the two indepen¬ 


dent variations of phases that preserve the minimum constraint (19). However, the step of the 
variation of each independent phase in coherent state counting must be larger than its quantum 
uncertainty. From the properties of the coherent states the phases satisfy the quantum uncer¬ 
tainty ^quant^k rsj 1/y/Nk- Thus, the step in variation of each phase cannot be smaller than this 
uncertainty. 


For the clarity of counting we should distinguish variation of the different phases. Notice that 


since the phases are constrained by the minimum condition (19), there are only two independent 
phase degrees of freedom, which can be chosen as 6q and 0 = — 6-±. The states of interest can 

thus be labeled as | ip) = \x,9o,0). 


Let us first fix the phase 6o and only vary the phase 6. The quantum uncertainly in this phase 
is A qua ntO ~ 1 / y/x and therefore for given x there are of order y/x independent coherent states 
that are orthogonal to each other. This set of states is obtained by variation of the relative phase 
6 with the elementary step bounded from below by the minimal quantum uncertainty W®- 
Now summing over different values of x over an interval 0 < x < x max , we get that, for fixed 9q, 
the total number of orthogonal coherent states crowded within the 1/N energy gap, is given by 
~ Xmax ~ y/~N ■ 


The story with the 8q phase is somewhat different. This phase is a Goldstone boson cor¬ 
responding to the generator Q o = diag( 1,—2,1), which is spontaneously broken already in the 
under-critical regime aN <C 1. Correspondingly, the coherent states that are orthogonal due to 
variation of 6q exist irrespective of criticality. The number of such orthogonal states for arbitrary 
value of x is ~ y/N — 2.x. Of course, they continue to contribute into the degeneracy also at the 
critical point, but this extra degeneracy has nothing to do with quantum criticality. 
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